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Abstract 


Stabilization is a fundamental requirement in the design of feedback compen- 
sators for flexible structures. The search for the largest neighborhood around 
a given design plant for which a single controller produces closed-loop stability 
can be formulated as an H^c control problem. The use of normalized coprime 
factor plant descriptions, in which the plant perturbations are defined as ad- 
ditive modifications to the coprime factors , leads to a closed-form expression 
for the maximum neighborhood boundary allowing optimal and suboptimal H 3 c 
compensators to be computed directly without the usual 7 iteration . This pa- 
per gives a summary of the theory on robust stabilization using normalized co- 
prime factor plant descriptions , and it describes the application of the theory to 
the computation of robustly stable compensators for the phase 0 version of the 
Control-Structures Interaction (CSI) Evolutionary Model . Results from the ap- 
plication indicate that the suboptimal version of the theory has the potential of 
providing the basis for the computation of low-authority compensators that are 
robustly stable to expected variations in design model parameters and additive 
unmodeled dynamics. 

Introduction 

In the design of controllers for physical systems, some trade-off is usually performed between 
design model accuracy and mathematical complexity. The more accurate analysis models often 
require computational time that is too excessive to qualify them as design models for control 
purposes. Also, many of the most widely used multivariable design techniques work best for 
moderate-order, linear, time-invariant design models (Maciejowski 1989). In practice, high-order 
nonlinear models are typically linearized about some operating condition and have their model 
order reduced to produce design models that conform to computational limitations or compensator 
implementation constraints. These practicalities introduce modeling errors in the form of unmodeled 
dynamics that must be accounted for in the controller design process. Additionally, parameters in 
the design and analysis models are not always accurately known and can cause destabilizing effects 
if parametric uncertainty is ignored or improperly treated. 

The foregoing considerations are especially critical in the design of controllers for flexible space 
structures (Joshi 1989). Space structure controller design models are generally obtained through 
some order-reduction procedure applied to a high-order analysis model obtained from finite element 
techniques. The order-reduction process essentially deletes a portion of the finite element model to 
produce a lower order controller design model. Although the unmodeled dynamics (represented by 
the deleted portion) are no longer contained in the design model, they can still be influenced by 
control inputs. Care must be taken in the design process to avoid control and observability spillover 
effects (Balas 1982) that destabilize the unmodeled dynamics. Also, the high-order model contains 
parametric uncertainties in natural frequencies, damping ratios, and mode shapes and slopes that 
get passed through to the design model. A fundamental requirement of control law design for flexible 
space structures is then the attainment and preservation of closed-loop stability in the presence of 
unmodeled dynamics and parameter uncertainties. 

Nonparametric uncertainties, such as unmodeled dynamics, typically occur in the high-frequency 
region of physical models. In controller design, unmodelcd dynamics are normally treated, often 
at the expense of reduced performance, by rolling off the compensator over a high-frequency band. 
Rolling off the compensator in this manner can be especially limiting for flexible structures in which 
unknown elastic modes are low-frequency, closely spaced, and fall within the bandwidth of the 
controller. The treatment of parametric uncertainties is still an open area of research (Dorato and 



Yedavalli 1990). In this study we are concerned with compensators that robustly stabilize a given 
system in the presence of both parametric and nonparametric uncertainties. 

Uncertainties may be viewed as perturbations about a nominal design model. If a single 
compensator stabilizes the nominal plant and, in addition, all systems within some neighborhood of 
the plant generated by the perturbations, the compensator is said to robustly stabilize the overall 
family of systems. The search for the largest neighborhood around a given design plant for which a 
single controller produces closed-loop stability can be formulated as an H x control problem. Glover 
and McFarlane (1989) show that the use of normalized coprime factor plant descriptions, in which the 
plant perturbations are defined as additive modifications to the coprime factors, leads to a closed- 
form expression for the maximum perturbation radius. The maximum radius can be computed 
directly in terms of the design model, thus allowing optimal and suboptimal robust compensators 
to be found without the usual 7 iteration of the H 00 design. 

This paper describes an application of the Glover-McFarlane robust stabilization theory to the 
control of a simulated structure configured to have many of the dynamic characteristics and controller 
design difficulties associated with flexible space structures. We begin with a description of the general 
robust stabilization problem followed by an overview of its results for normalized coprime factor 
plant descriptions. Mathematical background material is presented in appendix A, and the major 
computational algorithms for computing optimal and suboptimal robust controllers are presented 
in appendix B. Next, these theoretical results are applied to the computation of robustly stable 
compensators for the phase 0 version of the NASA Control-Structures Interaction (CSI) Evolutionary 
Model. Analysis of the compensators indicates that for the class of applications considered, the 
suboptimal version of the theory has the potential of providing the basis for the computation of 
low-authority controllers that are robustly stable to expected variations in design model parameters 
and additive unmodeled dynamics. 


General Robust Stabilization Problem 

In terms of transfer function matrix models of the design system, plant uncertainties (unmodeled 
dynamics and parameter inaccuracies) can be modeled in several ways. Let G(s) represent a p x m 
transfer matrix and AG(s) denote some perturbation to G(.s), both with real-rational elements. A 
perturbation is called additive if the perturbed plant (G A ) is written as 


and it is called multiplicative when 


or 


G a = G + AG 

(i) 

Ga = (Ip + AG)G 

(2) 

Ga = G(I m + AG) 

(3) 


A third method of modeling plant uncertainty involves the use of coprime factorizations (appen- 
dix A). Here, G is written in coprime factor form and the system perturbations are defined in terms 
of perturbations to the respective coprime factors. Any (stable or not) transfer function matrix 
(G(s)) can be represented in terms of a pair of asymptotically stable, real-rational, proper transfer 
function matrices that are coprime. For a left-coprime factor form with (M € £ RH » ), 


and the perturbed system (G A ) is given by 

G a = (M + AM )' 1 (N + AN) 


(5) 


2 



with 


( 6 ) 


AG = [AM, AN] € RHoo 


Controller design employing uncertainty models in equations (1) to (3) has been widely inves- 
tigated; see, for example, Chen and Desoer (1982) and Doyle and Stein (1981). Design methods 
employing coprime factor models such as equations (4) to (6) have not been as widely accepted 
or employed as the additive and multiplicative forms; however, they have been shown to have 
many theoretical and computational benefits by Vidyasagar (1985) and by McFarlane, Glover, and 
Vidyasagar (1990). 

Each uncertainty model can be represented as a special case of an upper linear fractional 
transformation 

Ga - ?u( P, AP) = P 22 + P 2 1 AP(I - P u AP)-‘ P 12 (7) 


where det(I — PnAP) / 0 and 


Pll 

P 12 

P 21 



! p 22 
1 J 


The standard plants P associated with three of the foregoing uncertainty descriptions are 


( 8 ) 


and 



AP = AG 


AP = AG 


O 

1 I 


t G 

! 

M -1 

! g 


l 


AP = [AN, -AM] 


(9) 


( 10 ) 


( 11 ) 


for equations (1), (2), and (5), respectively. Equation (7) represents a generalized uncertainty model, 
and the process of using feedback to stabilize and control P can be represented as the block diagram 
shown in figure 1. 


By employing figure 1, a general robust stabilization problem can be posed (McFarlane and 
Glover 1990). From viewing as a family of perturbed models for a given class of perturbations 
AP, one can seek a single compensator K (s) that stabilizes not only G (that is, G^ with AP = 0) 
but all members of the G^ family. If AP belongs to a class of admissible perturbations defined as 
the union of the set of stable bounded perturbations (RHoo) and the set of perturbations in RLqq for 
which G and G^ have an equal number of closed right half-plane poles, then the following theorem 
can be established (McFarlane and Glover 1990). 


Robust Stabilization Theorem 

For any P 22 of P given by equation (8) with stabilizable and detectable state- variable realization, 
the compensator K(s) of figure 1 stabilizes G^ = Fu (P,AP) for all admissible values of AP such 
that || AP Hoc < e if, and only if, 
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1 . K stabilizes G 


and 

2 . ||/|.(P,K) Ho, < e " 1 

where the lower linear fractional transform is given by 

T l { P, K) = P„ + P 12 K(I - P 22 K )- 1 P 21 (12) 

The parameter e in the theorem can be viewed as a measure of robust stability for a given closed- 
loop system. The problem of finding the largest level of robust stability is termed the optimal robust 
stabilization problem and is formally stated in the following discussion. 

Optimal Robust Stabilization Problem 

Find the largest strictly positive number e = e max such that for all admissible AP values satisfying 
|| AP Hoc < e, a single controller exists that stabilizes .F(/(P, AP). From the robust stabilization 
theorem, 

ei„ax= ||^(P,K)|| 0O ) (13) 

where K is chosen from all controllers that stabilize G. 

The computation of e max thus involves the solution of an H 0 c optimization problem such as that 
discussed by Francis (1987) or Doyle et al. (1989); that is, find 

^ ||^ ( P,K)|U=7,nm (14) 

over all controllers K that stabilize G. This H 0 0 problem can be posed for each representation 
of the additive, multiplicative, and coprime factorization uncertainty classes. In the additive and 
multiplicative cases, the solution typically involves a computationally intensive iterative procedure 
to find the smallest 7 (whereby 7 = 7 m in) such that the suboptimal robust stabilization problem 

|| J^(P,K) Hoc <7 (15) 

is solved. The solution of the optimal robust stabilization problem for the coprime factorization 
uncertainty class can also be approached in a similar manner; however, if the coprime factors of G 
are normalized, the 7 iteration procedure can be completely avoided and the computational effort 
greatly reduced. 

Robust Stabilization Problem for Descriptions of the Normalized Coprime 
Factor Plant 

Let G(s) have the coprime factor plant representation of equation (4). With P and AP given 
by equations ( 11 ), 

^£,(P,K)= * (I - GK )" 1 M ^ 1 (16) 

Also, because [AM, AN] e RHoc , the set of admissible perturbations contains only elements within 
RH^, whereby the need to observe the equality condition on the poles of G and G^ is completely 
eliminated. 
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When the coprime factors of equation (4) are normalized, that is, satisfy the relation 

M(s) M*(s) + N(s) N*(s) = I (17) 

the solution to the corresponding optimal or suboptimal robust stabilization problem has a surpris- 
ingly simplified form. For the problem of designing a controller K that robustly stabilizes a plant G 
written in normalized left-coprime factor form, the following properties hold: 

1 . A controller K is stabilizing for G and satisfies 

|| * (I - GK )' 1 M -1 ||oo < 7 = (! 8 ) 

if, and only if, K has a right-coprime factorization 

K = UV" 1 (19) 

for some (U € RHoo,V e RHoc) satisfying 

|| [— N, M] + y + ||oo < (1-7~ 2 ) 1/2 = (1-* 2 ) 1/2 ( 2 °) 

2. Solutions to the optimal robust stabilization problem using normalized left-coprime factoriza- 
tion give 

|| [*] (I-GK )- 1 M - 1 ||oc= {l- II [N,M] || ^| 1 (21) 

where the subscript H refers to the Hankel norm (which is discussed in appendix A). 

3. The maximum robust stability margin is 



4 . All optimal controllers are given by equation (19), where U and V satisfy 

|| [-N, M] + \yY ||o° = II [N,M] ||„ < 1 (23) 

Proofs of properties 1 to 4 may be found in either Glover and McFarlane (1989) or McFarlane and 
Glover (1990). 

Properties 2 and 3 state that for the robust stabilization problem using optimal normalized 
left-coprime factorizations (NLCF), a 7 iteration to solve the H problem is not necessary. The 
value of 7, n j n , given by the right-hand side of equation ( 21 ), can be foundthrough the computation 
of the largest Hankel singular value of the stable transfer matrix [N,M] obtained in the initial 
NLCF of G. Properties 1 and 4 show that both the optimal (eq. (14)) and suboptimal (eq. (15)) 
NLCF robust stabilization problems are solved by computing the closest (in the Hoc norm sense) 

[Ul* 

completely unstable (all poles in the open complex right half-plane) rational function y to the 
stable RHoc function |-N,M]. This computation involves a version of the Nehari extension problem 
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rui* 

(Francis 1987). By using y , U and V in RH^ can be constructed for use with equation (19) 

to form the compensator K. State- variable realizations for the subopt imal problem with 7 > 7 min 
can be found in Glover and McFarlane (1989) or McFarlane and Glover (1990). The theoretical 
approach to solve the optimal problem with 7 = 7 min can be found in Glover (1984). Algorithms 
for constructing state- variable realizations of the compensator K for both (optimal and suboptimal) 
problems are collected in appendix B of this paper. 


Note that the results obtained from application of the foregoing robust stabilization theory can 
be conservative. The approach of analyzing perturbations (AP) only in terms of their H 0 0 norm 
bounds rather than taking advantage of any known internal structure is fundamentally conservative 
by nature (Maciejowski 1989). Known structural properties of perturbations may be incorporated 
by using structured singular value techniques such as those implemented by Balas et al. (1991). In 
addition, a feedback controller K may well stabilize G and either for some particular structured 
AP such that || AP ||oo > e m ax or for all structured AP terms within a subset of || AP ||* < e max 
containing AP = 0 while not satisfying the more restrictive condition of simultaneously stabilizing 
all elements of G^ for || AP ||* < e max . 


Loop-Shaping Procedure Within the NLCF Robust Stabilization Structure 

Because G = M -1 N is a normalized left-coprime factorization, we can write 


(I-GK)-'M 




K 

I 


(I - GK) -1 M -1 [M,N] Hoc 


K(l-GK)" 1 
(I - GK) 1 


K(I - GK)-‘G 

(I — GK) _1 G 


< 


7 


(24) 


Thus, for either the optimal or suboptimal versions of the NLCF robust stabilization problem, the 
satisfaction of equation (18) for 7 > 7 min also implies that (Safonov and Chiang 1988) 


and 


HK(I-GK)- 1 Hoo < 7 

(25) 

|| (I — GK) -1 H* < 7 

(26) 

||K(I-GK) -1 G ||oo < 7 

(27) 

|| (I-GK)-'G ||* < 7 

(28) 


The significance of the individual transfer function matrices used in equations (25) to (28) may be 
seen from analyzing the feedback control loop in figure 2. The system G is subject to exogenous 
inputs in the form of reference commands or actuator noise (r), disturbances reflected to the output 
(d), and measurement noise (7/). From figure 2 we can derive 


y( s ) = (I - GK)->GK r,(s) + (I - GK)“ l G r(s) + (I - GK)- 1 d(a) (29) 

u(s) = K(I - GK) -1 [d(s) -fi r](s)] + (I - KG) -1 r(s) (30) 

u(s) - r(s) = K(I - GK)- 1 |d(s) + i|(a)] + K(I - GK) ‘G r(s) (31) 

In equations (29) to (31) the standard notation in writing f(s) to represent the Laplace transforma- 

tion of /(() has been used. Equations (29) to (30) reflect the well-known inherent trade-off between 
attenuation of disturbances (d) and filtering out measurement and/or actuator noise ( 77 , f). Be- 
cause the sensitivity ((I - GK) -1 ) and complementary sensitivity ((I - GK) _1 (-GK)j transfer 
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function matrices are related by 

(I - GK)" 1 + (I - GK) -1 (-GK) = I (32) 

both cannot be kept small over the same frequency range. Normally, the disturbances are assumed to 
be large in magnitude only over a low-frequency range, and actuator and sensor noise are appreciable 
only over a complementary high-frequency range. The compensator K is typically designed to cause 
|| (I - GK) 1 ||oo to be small over low-frequency ranges and || (I - GK) _1 GK Hoc to be small over 
high-frequency ranges. 

The NLCF robust stabilization compensators cause equations (25) to (28) to be satisfied over 
the entire frequency range, whereby the output y may be subject to any high-frequency noise. Such 
design difficulties are normally addressed through frequency-dependent weighting matrices applied 
directly to the sensitivity and complementary sensitivity matrices before compensator design. Only 
a special class of weighting matrices are shown to be allowed if the exact-solution advantage of the 
NLCF robust stabilization problem is to be preserved (McFarlane and Glover 1990). 

Let W,(») and W 2 (.) be system precompensator and postcompensator matrices, respectively, 
and define an augmented plant G A (s) by 


g a( s ) = W 2 (s) GO) W\{s) (33) 

Performing an NLCF robust stabilization design with G replaced by G^ yields a dynamic com- 
pensator K A (s) that robustly stabilizes G^. Figure 3(a) gives a block diagram illustrating this 
loop-shaping procedure. Simple block manipulation yields figure 3(b), which shows that the corre- 
sponding compensator K to be applied to the unshaped plant G is 

K(s) — Wj 0) K a (s) W 2 (s) (34) 

We then have 


l, max = II [N A ,M A ] I&) 


1/2 


inf 

K 


inf 

K A 


Wf 1 K(I-GK)- 1 W 1 
W 2 (I — GK) -1 Wj 1 


I 


(I-G^)- 1 !!^] 


Wf'Ktl-GK^GW! 
W 2 (I - GK) _1 GWj 


(35) 


which indicates the weighting configuration that can be applied to the elements of equation (24) 
if the exact-solution structure is to be preserved. In general, if other weighting configurations are 
desired, the normal 7 iteration procedure is required. 

The introduction of dynamic weighting matrices always increases the order of the compensator. 
Let n be the order of a minimal realization of G(s), m be the order of W^s), and n 2 be the order 
of W 2 (s). Given state-variable realizations for Wi(s), W 2 (s), and K^(s), repeated application of 
equation (A5) yields a state-variable realization for K(s). The order of the realization for K(s) will 
be < n + 2 {n\ -I- n 2 ) after elimination of any uncontrollable and unobservable modes. Lower order 
compensators can be obtained through application of compensator order-reduction methods such as 
that described by Anderson and Liu (1989). However, compensator order-reduction methods should 
be used with caution because they can attenuate robust stability. 


Stability Robustness to Unmodeled Dynamics 

In flexible space structures, perturbations in the form of unmodeled dynamics occur naturally in 
the additive form (eq. (1)). A direct relationship does not seem to exist in defining [AN, AM] as a 
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function of a given additive AG in equation (1). Setting 

G a = M -1 N + AG = (M 4- AM) -1 (N + AN) 

yields 


(M + AM)AG = [AN, AM] 
For small || AMAG ||oo, equation (37) becomes 

AG = M _1 [AN, AM] 


Im 

-G 


Im 

-G 


(36) 

(37) 

(38) 


for which a solution is 

[AN, AM] = M AG(I m + G*G) _1 (I m , -G*) (39) 

Then, robust stability is preserved if 

|| M AG(Im + G G) l(I m ,G ) ||oo < e max 

Rigorous use of equation (39) is likely to produce conservative designs. In the application of the 
NLCF robust stabilization theory that follows, stability robustness to additive unmodeled dynamics 
in the form of equation (1) is incorporated by use of the weighting matrices in equation (34) and is 
analyzed through the examination of closed- loop conditions for stability robustness (Francis 1987) 

|| AG K(I-GK)- 1 ||oc < 1 ( 4 °) 

or the more conservative sufficient condition for inequality (40) 

ct[AG(s)]ct[K(I - GK) _1 (s)] < 1 (41) 


for s = joj and u> £ [0, oo). 


In figure 2, let 


and 


K(s) = 


G(s) = 


A c 

B c 

Cc 

D c 

A 

B 

C 

D 


(42) 


(43) 


Then, state- variable realizations for the transfer matrices in equations (29) to (31) can be formed, 
wherein each has the same system matrix given by 


1” A + B(I — D c D) _1 D c C B(I — D c D) -1 C c 1 
Ac ' = [ B C (I — DD c ) _1 C A c -I- B C D(I — D c D) _1 C c J 

An eigenvalue analysis of A ^ often gives an indication of the effect of the compensator on the open- 
loop eigenvalues of A. An eigenvalue analysis of A c / can occasionally give a measure of stability 
robustness to additive perturbations when the realization (A, B, C, D) in equation (44) is replaced 
by a state- variable realization of G^. 
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Robust Stabilization of the Phase 0 Evolutionary Model 

In this section , the previously discussed theory and algorithms for the robust stabilization of a 
plant modeled in normalized left-coprime factorization form are applied to produce compensators 
for the control of a model of a laboratory structure that has many of the characteristics and design 
difficulties associated with flexible space structures. Except for certain high-order transient response 
simulations, all computations were performed using a 486 personal computer employing the 386 
MATLAB computational environment in Anon. (1990). Use was also made of software from the 
MATLAB Control System Toolbox (Little and Laub 1986) and Robust-Control Toolbox (Chiang 
and Safonov 1988). Balanced realizations were computed by using the algorithm based on Singular 
Value Decomposition described by Laub et al. (1987). The H ^ norms were calculated by using the 
bisection algorithm of Boyd, Balakrishnan, and Kabamba (1989). 

Description of Phase 0 Model 

The Control-Structures Interaction (CSI) Evolutionary Model is a laboratory testbed concept in 
which a sequence of testbeds is evolved with each new facility having more challenging dynamics 
and control characteristics than the previous one. The testbeds are to be designed and constructed 
at the Langley Research Center for the experimental validation of control techniques and integrated 
design methodology developed under the CSI program at Langley (Newsom et al. 1990). The phase 0 
model was the first testbed to be constructed under this program, but, unfortunately, the phase 0 
model is no longer in existence at Langley. However, many useful studies were performed using the 
phase 0 model (for example, see Lim, Maghami, and Joshi (1992); Maghami, Joshi, and Armstrong 
(1993); Lim and Balas (1992), and its data base has been archived and is still available for this and 
future studies. A schematic of the phase 0 structure is shown in figure 4, and a detailed description 
can be found in Belvin et al. (1991). 

The phase 0 structure consisted of two vertical towers and two horizontal booms attached to a 
central 62-bay-truss main section with each ,bay being a 10-in. cube. (See fig. 4.) The structure 
was suspended from the laboratory ceiling by springs and two long cables designed to minimize the 
interaction between the suspension and the structural modes. A laser source was mounted at the top 
of the taller vertical tower, and a 16-ft reflector with a mirrored surface was mounted on the shorter 
tower. The laser beam was reflected by the mirrored surface onto a detector surface 660 in. above the 
reflector. The total structural weight was 741 lb. Eight proportional bidirectional gas thrusters (air 
jets) provided the input force actuation, and eight nearly collocated servo accelerometers provided 
output measurements. 

Global line-of-sight pointing studies using the laser targeting system have been performed by Lim 
and Balas (1992). The present study is directly concerned with only vibration suppression about 
a given operating point and does not treat laser targeting as such. However, vibration suppression 
of the laser tower modes is a critical part of our design requirements because damping of the laser 
tower structure was a troublesome issue in previous laboratory tests. 

The phase 0 model had six nonstructural modes (due to suspension) and many significant elastic 
modes. The NASTRAN finite element model (Belvin et al. 1991) consisted of 3560 degrees of 
freedom. A total of 86 modes with frequencies below 50 Hz were selected as a truth or evaluation 
model. For this study, as in previous studies (Lim, Maghami, and Joshi 1992), a reduced-order 
model consisting of 25 modes (selected from the 86 modes through a controllability and observability 
analysis) was used for the controller design model. Table I shows the frequency range of these 
25 modes in hertz. Frequencies of the first six pendulum/suspension modes, brought about by the 
cable suspension in a 1 g environment, range from 0.147 to 0.874 Hz. The first two elastic modes 
(7 and 8) are truss beam bending in the XY- and XZ-planes (defined in fig. 4) with frequencies of 
1.474 Hz and 1.738 Hz, respectively. A uniform damping ratio of 0.5 percent is assumed for all the 
modes. 
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By using data from the finite element analysis, a dynamic mathematical model in the modal 
coordinate system can be constructed. A 50th-order state-variable realization of this model in 
(block) phase-canonical form will appear as (Ay, By, Cy,Dy) where 


A y = 





Cy = [0,$]Ay Dy = ** T 


(45) 


mtii / a s?\ 

fl = diag(uq , ui 2, ^25) ( 46 ) 

Z = diag(2(iuq, 2C2^2, 2^25^25) 

and $ is an 8 x 25 matrix of mode shapes obtained from a finite element analysis (Belvin et al. 
1991). In equations (46) and (47), u>j denotes frequency and Ct denotes a damping ratio of 0.005 
for i = 1, ... ,25. Eigenvalues of Ay are given in table I. Because the damping ratios are small, the 
frequencies in radians per second are closely approximated by the imaginary parts of the eigenvalues. 
Figure 5 shows an open-loop frequency response plot of the singular value bounds of the transfer 
matrix Gy defined by the (Ay,By,Cy,Dy) system. 

In flexible structures, higher frequency modes are more difficult to measure and compute 
accurately. For the phase 0 structure, the finite element model provided reasonably accurate natural 
frequency and damping values for modes below 2 Hz. However, higher frequency modes, beginning 
with the 10th mode at 2.301 Hz, are not accurately known (Lim and Balas 1992). In this design, 
only modes with frequencies up to 14 rad/sec are used to form the compensator design model, and 
modes with higher frequencies are used to represent unmodeled dynamics. Therefore, modes above 
mode 9 are truncated from the 25-mode model and are accounted for as an additive uncertainty in 
the design process. The matrix Gy now appears as in equation (1) with 

Ga = Gy = G + AG (48) 

G = C(sIi8-Ar 1 B + D (49) 


AG = Cf(sl 3 2 — A() _1 Bf -I- D< (50) 

The realizations (A,B,C,D) and (A f , B t , C t , D t ) are obtained by collecting and rearranging 
appropriate rows and columns of (Ay, By, Cy, Dy). Numerical data for G and AG are given in 
appendix C, and figure 6 shows individual frequency response plots for G and AG. 


Design Objectives 

The objective of the control system design is to increase the damping of all the pseudo rigid-body 
and structural modes of the design model G. The designs must also possess stability robustness with 
respect to unmodeled structural modes, of which AG from equation (50) is taken as a representative 
sample, and, if possible, possess stability robustness to parametric uncertainties such as errors 
expected in frequency, damping, and mode shapes in the design model. 

Loop-Shaping Procedures 

By using the numerical data for the 18th-order system (A, B, C, D) found in appendixC, the 
construction given by equations (A18) to (A22) leads to a state- variable realization for [N,M] with 
a Hurwitz system matrix. Computing the Hankel norm of [N,M], as outlined in appendix A, gives 
|| [N,M] ||h = 0.8972, whereby, from equation (22), e ma x = 0.4417 and 7min = l/e m ax = 2.264. 
An optimal compensator (using eq. (B16) with r = 1) and a central suboptimal compensator (using 
eq. (B28) with e = 0.9e max ) were obtained following the algorithms in appendix B. Frequency 
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response plots for the compensators are shown in figure 7. Both compensators, when applied to G 
in the feedback fashion of figure 2, enhanced the stability of G but grossly violated the condition 
stated in inequality (40) with AG given by equation (50). Violation of inequality (40) indicates that 
the compensator 

u (s) = K(s) y(s) (51) 

applied to G^, given by equation (48), has the potential of a destabilizing spillover effect (Joshi 
1989) on the AG dynamics. The failure to satisfy the condition in inequality (40) was primarily 
caused by the lack of free parameters for adjustment in the algorithms and the fact that G(s) is not 
strictly proper, in which case both optimal and suboptimal compensators will not be strictly proper 
as can be seen from figure 7. 

The spillover problem could be resolved through incorporation of loop-shaping functions of the 
form given in equation (33). The weighting functions employed were 


WiW = i 8 

(52) 

and 


W 2(5) = W(s) - f(s) I 8 

(53) 

where the scalar function f(s) is given by 


/(s) = k/(s + a) 1 

(54) 


The positive real parameters a,z, and k are adjusted from observation of the condition given by 
inequality (41) with the compensator in equation (34). The parameter a is chosen such that a plot 
of the inverse of a[K(I - GK) -1 (ju>)] versus u> € [0, oo) breaks upward before ll> = 14 rad/sec, the 
approximate frequency at which the AG dynamics become predominant. The parameter i roughly 
controls the slope of the upward break and was taken as i — 1 or 2. The quantity k adjusts the 
magnitude of W(s) and ranges between 0.08 and 4.0 in the following discussion. Minimal-order 
state- variable realizations for W(s) for i = 1 are 



where 


and, for i — 2, are 


where 



0 0 ... 0 0 
0 0 ... 0 1 

1 0 ... 0 0 

0 0 ... 0 0 




(56) 


(57) 


( 58 ) 
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(59) 


and 



c[? = (o 8x8 | *b$) 


(60) 


D$ = 0 8x8 


(61) 


Realizations of equations (55) to (61) were obtained from the Control System Toolbox function 
entitled u tf2ss.m. ,? 


Increasing the parameter k not only moves the closed-loop eigenvalues of the controllable modes 
of the design model farther into the complex left-hand plane but also increases the potential of 
destabilizing spillover into the modes of A*. Observing the real parts of modes 10 and 11 shows that 
these modes need stability augmentation as much as modes 1 to 9. However, of the first 11 modes, 
modes 10 and 11 are the least controllable (with mode 11 more controllable than mode 10). A 
Hankel singular value analysis of the 11-mode model yields the mode numbers ordered in decreasing 
controllability and/or observability as 8, 7, 3, 2/9 (tie), 6, 4, 1, 5, 11, and 10. Including modes 
10 and 11 in the design model and increasing k to provide the added stability invariably produces 
compensators (optimal and suboptimal) that violate conditions given in inequalites (40) and (41). 
By using the methodology and weighting described herein, modes 10 and 11 will necessarily remain 
relatively unchanged when stability to additive perturbations is a design requirement. 


Robustness to Structured Perturbations 

For modes 1 to 9 below 2 Hz, the accuracy of values of natural frequency and damping ratio are 
within 1 percent and 10 percent, respectively. (See Lim and Balas 1992.) In order to evaluate 
the compensators for perturbations in frequencies and damping ratios within these ranges, the 
frequencies and damping ratios Q in the A matrix of the 9-mode design model were replaced 
by perturbed values (c j p )i and (Cp)i given, respectively, by 


( u >p)i - Wj -f (Suj)iUJi (62) 

and 

(Cp)i = Ci + (SOiQ (63) 

for i = 1,...,9. In equations (62) and (63), ( 6u;)i and (6£)i are random variables uniformly 
distributed within [-0.01, 0.01] and [—0.1, 0.1], respectively. The new perturbed system matrix 
is denoted by A$. If G (s) is given by equation (49), the transfer matrix for the perturbed system is 
given by 

g 6( 5 ) = C*(*Iig - A^) -1 !* + D (64) 

where 

Cfi = CA-'Ai (65) 

to reflect acceleration measurements. 

With 

G a ( S ) = W 2 (s) G(s) Wi(*) 

and 

G 6 a (s) = W 2 (s) G 6 (s) Wj(s) 
normalized left-coprime factors are found such that 

Ga(s) = (M^)" 1 ^ (68) 


( 66 ) 

(67) 
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and 


(69) 


G«» = 

Define __ 

AP = [AN,,,-AM„] (70) 

where ___ 

AM a = M 6 a -M a (71) 

and 

AN4 = - N.4 (72) 

If || AP Hoc < e A < M.maxi the compensators also stabilize the perturbed system in addition to 
stabilizing G(s). 

Compensator Design 

Optimal compensator . The parameters (A:, i, a) in equation (54) were adjusted to meet the 
design objectives. Because the optimal compensator was not strictly proper, i = 2 was found to 
provide faster roll-off and to best allow the satisfaction of conditions in inequalities (40) and (41). 
A representative optimal compensator used (fc,i,a) = (0.1, 2, 0.5). After scaling each channel 
of G by 0.1 /(a + 0.5) 2 , the result was e^max = 0.6749. Figure 8 shows the scaled version 
of G(s), and figures 9 and 10 show singular value bounds for the corresponding compensators 
K^(s) and K(s), respectively. Table II gives the eigenvalues of the Hurwitz system matrix from 
a state- variable realization of K^(s). The order of the compensator K A is 33 (corresponding to 
n = 18, m = 0, n2 = 16, r = 1, and i = 2), whereby the final compensator K is on the order of 49. 

Over 2000 realizations of equations (62) and (63) were computed for the optimal compensator and 
corresponding AP transfer matrices (given by eq. (70)) tested for satisfaction of || AP || oo < max- 
No violations were encountered. Values of || AP ||oo ranged between 0.1728 and 0.5888 with a mean 
of 0.4627, and the standard deviation was 0.0726. Satisfaction of conditions in inequalities (40) 
and (41) is depicted in figure 11, which indicates that mode 20 (the laser tower mode) at 41.9 rad/sec 
(6.7 Hz) is the mode most likely to experience destabilizing spillover. (In this report, frequency is 
expressed in radians per second or Hertz (or in both) for comparison purposes with the literature.) 
This property has also been observed experimentally in previous studies. The peak value f ffx norm) 
of the curve in figure 11(b) is 0.024 at 41.9 rad/sec, which indicates an additive stability robustness 
margin of about 97 percent. 

This ultraconservative margin for additive stability robustness was forced by the desire to also 
have a compensator that guarantees stability robustness to expected parametric uncertainties in 
frequencies and damping ratios. If k is increased to 2.0 for the same values of i and a, the peak 
value of the curve in the corresponding figure 11(b) is 0.570, which indicates a more reasonable 
margin of 43 percent. However, k = 2.0 produced e Amax = 0.4856, and this reduced value of e Amax 
leads to the occurrence of violations greater than 40 percent of the time in the random tests for 
|| AP Hoc < e^ max . For fixed values of i and a, decreasing k increases e A>max and decreases 
|| AP || oo. A value of k allowing a sufficiently wide “gap” between the two quantities is required 
when specific ranges of parametric variations are considered. 

A state- variable realization ( A c , B c , C c , D c ) of the optimal compensator K was formed, and the 
controller was applied (in the manner of fig. 2) to the 50th-order system G^. The closed-loop system 
matrix appears as A c / in equation (44), but with (A, B, C, D) replaced by (Ay, By, Cy,Dy). 
When comparing the imaginary parts of the eigenvalues of A c / with the imaginary parts of the 
open- loop eigenvalues (given in table I), certain A c / eigenvalues can be identified as the closed-loop 
modifications of the eigenvalues of Ay. These eigenvalues and corresponding damping ratios (which 
are correlated with open-loop mode numbers) are given in table III. The compensator leaves the 
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open-loop eigenvalues essentially unchanged except for modes 1 to 3. Optimal compensators with 
better stability augmentation can be obtained at the expense of violations in robust stability because 
of parametric variations in frequencies and damping ratios of the design model. 

Suboptimal compensators . Suboptimal compensators computed from strictly proper systems 
G>t(s) are strictly proper, whereby lower order weighting functions (other than those used in 
the foregoing optimal compensator studies) can be employed without difficulty in satisfying the 
additive robustness conditions. Suboptimal compensator studies were performed using i = 1 and 

= 0.9e4 max from which a representative result had (fc,i,a) — (0.5, 1,0.1). After scaling each 
channel of G by 0.5/(s + 0.1), the result was found to be = 0.667. Figure 12 shows 

the scaled version of G, and figures 13 and 14 show singular value bounds for the compensators 
K^(s) and K(s), respectively. Table IV gives the eigenvalues of the Hurwitz system matrix from 
a state variable realization of K^(s). The order of the compensator is 26 (corresponding to 
n = 18, n\ — 0, ri2 — 8, and i — 1), whereby the order of the final compensator K is 34. 

Again, over 2000 realizations of equations (62) and (63) were computed for the suboptimal 
compensator and the corresponding AP transfer matrices tested for satisfaction of || AP ||oo < 
e A, max no violations encountered. Values of || AP ||oo ranged between 0.0561 and 0.5595 with 

a mean of 0.2798 and a standard deviation of 0.0907. The expected variations in frequency and 
damping ratio of 1 percent and 10 percent, respectively, are apparently close to the upper bounds 
for robust stability for k = 0.5. Increasing these variations to 1.5 percent and 15 percent causes 
violations about 6 percent of the time. 

Satisfaction of conditions in inequalities (40) and (41) is shown in figure 15. An interesting note 
is that figure 15(a) indicates that mode 20 (the same mode as that indicated by fig. 11) is the mode 
most likely to be troublesome, whereas figure 15(b) flags mode 15 at 27.6 rad/sec. The peak value 
of the curve in figure 15(b) is 0.562, which indicates an additive stability robustness margin of about 
44 percent. 

The suboptimal compensator was also applied to the control of the 50th-order system G^, and 
an eigenvalue analysis was performed on the resulting matrix. The results are given in table V. 
The fact that the real parts of the eigenvalues for modes 10 to 25 are not significantly changed from 
the corresponding values of table I indicates that the additive robustness conditions are satisfied. 
Modes 3 to 9 all have enhanced stability except for possibly mode 5. The eigenvalue data give two 
eigenvalues with imaginary parts close to the imaginary part of open- loop mode 5. Both are shown 
in table V. The entry with the largest real part is possibly the closed-loop eigenvalue of mode 5 
because this mode is the least controllable and observable of the first nine modes. No correlation 
could be made for modes 1 and 2. 

For comparative purposes, a calculation was performed with (fc,z,a) = (0.5, 1,0.1) and e ^ = 0. 
Compensator equations for this calculation are given by the 7 — - ► 00 case found at the end of 
appendix B. Figures 16 to 18 and tables VI and VII provide the same information as previously 
provided for the optimal and other suboptimal compensators. Figure 18 indicates that robustness to 
additive perturbations measured by conditions in inequalities (40) and (41) is somewhat improved. 
The peak value of the curve in figure 18(b) is 0.184 at 27.6 rad/sec with a robustness stability margin 
of about 81 percent. 

For fixed values of the parameters as increased from 0 toward max within the 

suboptimal structure, the magnitudes of the imaginary parts of the eigenvalues of a state-variable 
realization of K^(s) generally decreased, whereas the magnitudes of the real parts generally 
increased. Thus, for small values of a relatively slower compensator was produced; this caused 
smaller values on the left side of inequality (41) which, in turn, created larger values of additive 
robustness stability margin. 
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Evaluation model simulations . The four compensators discussed were applied to the control 
of the full 86-mode-evaluation phase 0 model subjected to a transient input disturbance. Although 
not directly considered in the robust stability design process, this application gave an indication 
of the relative merit of the compensators in the area of closed-loop performance. The simulation 
consisted of (1) applying an excitation input sequence for the first 9 sec, (2) allowing 1 sec of free 
decay, and then (3) applying a controller at the 10-sec mark for a total duration of 30 sec. The 
input sequence consisted of harmonic forces designed to excite two pendulum modes (1 and 3) and 
the first two bending modes (7 and 8) using a single actuator for each mode. Specifically, actuators 
1 (mode 7), 2 (mode 8), 6 (mode 3), and 7 (mode 1) were excited with signals of 1.474 Hz, 1.738 Hz, 
0.155 Hz, and 0.147 Hz, respectively. No actuator dynamics were considered. 

The compensators were discretized at a sampling rate of 133 Hz. The input sequence and 
sampling rate are the same as those employed by Lim, Maghami, and Joshi (1992) in their 
experimental investigation using the actual phase 0 structure. The format for the presentation 
of the results herein parallels that of the experimental study. Figures 19 to 34 show input and 
output responses for collocated thrusters and accelerometers 7 and 8 (laser tower) for the optimal 
and suboptimal controllers. The best performance is obtained by the suboptimal compensator with 
M.max = 0.667 followed closely by the suboptimal compensator with = 0. The optimal (2.0, 2, 
0.5) compensator had marginal performance characteristics, and the highly conservative robust- (0.1, 
2, 0.5) compensator had almost no effect. The suboptimal responses compare favorably with those 
of Lim et al. which were designed with performance issues in mind. 

Concluding Remarks 

In this report the application of a robust stabilization methodology to the control of a simulated 
model of an actual laboratory apparatus has been considered. The apparatus, known as the Control- 
Structures Interaction (CSI) Phase 0 Evolutionary Model, was built at the Langley Research Center 
for the purpose of investigating some of the problems and solution approaches involved in the control 
of flexible space structures. 

With the view that system uncertainties are perturbations about a nominal design model, a single 
robustly stabilizing compensator was sought that not only enhanced the stability of the nominal 
plant but, in addition, stabilized all systems within some neighborhood of the plant generated 
by the perturbations. The computation of the largest neighborhood around a given design plant 
for which a single controller produces closed-loop stability was shown to be formulated as an 
Hoc control problem. Several formulations of the H 0 c problem were possible depending on the 
character of the perturbation space. A formulation based on the G lover- McFarlane theory (1989) 
was employed wherein descriptions of the normalized coprime factor plant were used and the plant 
perturbations were defined as additive modifications to the coprime factors. An attractive feature 
of the normalized coprime factorization approach is that it leads to a closed-form expression for 
the maximum perturbation radius. The maximum radius can be directly computed in terms of the 
open-loop design model, thus allowing optimal and suboptimal robust compensators to be found 
without the usual 7 iteration of the H x design. 

Computational algorithms for both the optimal and suboptimal versions of the McFarlane-Glover 
theory (1990) were summarized, and representative compensators were computed and analyzed 
for the control of the simulated Langley phase 0 structure. Through incorporation of weighting 
functions, the compensators were made to roll off in such a way as to cause the closed-loop system 
to be robust to additive unmodeled dynamics. A Monte Carlo approach was employed to test the 
robust stability of the compensators against expected structured parametric perturbations in the 
compensator design model. Even though these calculations offer no mathematical proof, they do 
give a degree of confidence that the compensators robustly stabilize the nominal system for expected 
perturbations in frequencies and damping ratios of the system matrix of the design model. The 
compensators were tested through application to a full-order evaluation model. 
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Results from the study indicate that when requiring the compensators to satisfy all design ob- 
jectives of stability augmentation, robust stability to unmodeled dynamics appearing as additive 
perturbations, and robustness to structured parametric variations, the optimal robust compensators 
can be overly conservative with marginal stability augmentation, whereas the suboptimal compen- 
sators are not. For the class of flexible structure applications considered, the suboptimal version 
of the McFarlane-Glover theory provides a viable approach for the computation of low-authority 
controllers providing robust stability augmentation for variations in design model parameters and 
additive unmodeled dynamics. These controllers may need to be supplemented with high-authority 
loops to provide additional performance. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
March 29, 1993 
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Appendix A 


Nomenclature and Definitions 

This appendix offers a compilation of certain fundamental results from linear systems theory that 
are employed in the body of this report. 

State- Variable and Transfer Matrix Representations 

Assume that a dynamic system is modeled by a real, linear, finite-dimensional, time- invariant, 
continuous set of ordinary differential equations written in vector-matrix form as 


x(t) = A x(t) + B u(t) 
y(0 = C x(t) + D u (t) 


where, for t G [0,oo), 


x(0 G R n 
y(t) € R p 
u(f) G R m 
A G R nxn 
B G R nxm 
C G R I>X " 
D G Ri ,xm 


state vector 
output vector 
control vector 
system matrix 
control-effectiveness matrix 
output matrix 
control-feedthrough matrix 


(Al) 

(A2) 


The transfer /; x rn function matrix (G(,s)) for the (A, B, C, D) system is given for ,s G C by 


G(.s) = C(.sl - A) _1 B + D : = 
If D is square and nonsingular, 

g-‘oo = 


A 

B 

c 

D 


A - BD ’C 

BD 1 

-D -1 C 

D 1 


if 


G, = 


Go - 


Ai 

CQ 

Ci 

D, 


A 2 

b 2 

c 2 

d 2 


(A3) 


(A4) 


then state- variable realizations (not necessarily minimal) for the cascaded system G 1 G 2 are 
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Ai B1C2 B1D2 

G1G2 - 0 a 2 b 2 

Cl D!C 2 D1D2 


(AS) 


where 0 represents a null matrix of appropriate size. 

When A is Hurwitz, that is, has eigenvalues only in the open, complex left half-plane, the 
reachability gramian (P) and the observability gramian (Q) are given (Kwakernaak and Sivan 1972) 
as the unique symmetric nonnegative definite solutions to the Lyapunov equations 

AP + PA r + BB t = 0 (A6) 

A r Q + QA + C T C = 0 (A7) 

If (A, B) and (C, A) are, respectively, completely controllable and completely observable, then P 
and Q are symmetric and positive definite. 

Denote the ith eigenvalue of a square matrix E by A Z (E). When Re[A 2 (A)] < 0 for i = 1, 2, . . . , n, 
the Hankel singular values of the transfer function matrix G(s) are defined as (Glover 1984) 

*i[G(a)] = {Ai(PQ)} 1 / 2 (A8) 

where P and Q satisfy equations (A6) and (A7), respectively, and, by convention, 

°i[G(s)] > <r i+1 [G(s)] (A9) 

for i - 1,2, . . . , n — 1. The Hankel norm of G(s), denoted by || G ||//, is defined as 

II G || h = ^i[G(s)] (A10) 

which is the largest Hankel singular value of G(s). The eigenvalues of the product PQ are invariant 
under coordinate transformations on x(f) in equations (Al) and (A2). When A is Hurwitz, a 
particular coordinate transformation called a balancing transformation exists (Moore 1981) such 
that 

P = Q = diag(<7 1 , rr 2 , . . . , cr n ) (All) 

Computational aspects of balancing transformations are discussed in Laub et al. (1987), and 
applications to model-order reduction are surveyed in Anderson and Liu (1989). 

Let RLr x; denote the space of all real-rational proper transfer function matrices that have no 
poles on the imaginary (s = juj) axis of the complex plane and RH^ C RLr ^ denote the subset of 
all asymptotically stable (no poles in the closed right half-plane) proper transfer function matrices. 
With G(s) 6 RHoc, the //oc norm of G(,s) is 

II G Hoc = H " P ' v[G(jv)] (A12) 

u) E R 

where a denotes the largest singular value (not the Hankel singular value) of the complex matrix 
G (ju;) for a given real uj. For a discussion of matrix singular values and their applications in linear 
systems theory, see Golub and Van Loan (1989) and Kleiria and Laub (1980). Computational aspects 



18 



of Hoc norms are discussed in Boyd, Balakrishnan, and Kabamba (1989) and Boyd and Balakrishnan 
(1990). 

Coprime Factor Representations 

In addition to state- variable and transfer-function matrix representations, all time-invariant linear 
systems have coprime factor representations (Vidyasagar 1985). Suppose that M G RH oc an d 
N G RH oc have the same number of rows. Then M and N are left coprime if, and only if, Ur RH oc 
and WeRH oc exist such that 

M (s)V(s) + N (s) U(s) = I p (A 13) 

for all s G C. Equivalently, M and N are left coprime if, and only if, [N, M] has a right inverse 
in RH oc> Any (stable or not) transfer function matrix G(s) can be represented in terms of a pair 
of asymptotically stable, real-rational, proper transfer function matrices that are left coprime. This 
representation, termed a left-coprime factorization of G(s), is given by 

G(s) = M(s) _1 N(s) (A14) 

where M G RH 0 c is square with det(M) / 0, and N G RH 0 c and M are left coprime. A particular 
left-coprime factorization, called a normalized left-coprime factorization (NLCF), is one in which 


N(s)N*(s) + M (s) M » = Ip (A15) 

for all s G C where 

N*(.s) = N r (-s) ( A16) 

M*(s) = M^(— s) (A17) 

A state- variable realization for a normalized left-coprime factorization can be formed from a minimal 
(A, B, C, D) realization of G(s) (Vidyasagar 1988). First, solve the generalized filter algebraic 
Riccati equation (GFARE) 

(A - BS *D r C)Z + Z(A - BS _ 1 D r C) r - ZC r R _1 CZ + BS 'B 7 ' = 0 (A18) 


for the unique symmetric positive definite matrix Z where 

S = I w + D r D (A19) 


Then, with 


and 


R = I p + DD 7 

H = ~{ZC T + BD 7 )R 1 


[N | M] 


A + HC 

B + HD 

H 

pa 

i 

To 

o 

r-'/ 2 d 

R -’/‘ 2 


( A20) 

(A21) 


(A22) 


N and M form a normalized left-coprime factorization of G(.s) constructed as in equation (A 14). 

Similar to left-coprime factorizations, all rational transfer matrices also have right -coprime 
factorizations (Vidyasagar 1985). Suppose that M G RH 0 c and N G RH oc have the same number of 
columns. Then, M and N are defined to be right coprime if. and only if, U G RH o c and V G RH oc 
exist such that 

VWM(s)+U(s)NW = I fH (A 23) 
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for all s E C. That is, the matrix 


' N * 

M 


has a left inverse in RHoq. The transfer matrix G(s) is said 


to have a right-coprime factorization when 


G(s) = NM" 1 


(A24) 


where M € RHoq is square with det(M) ^ 0, and N e RH ^ and M are right coprime. A particular 
right-coprime factorization, called a normalized right-coprime factorization (NRCF), is one in which 

M*(s) M(s) -f N *{s) N(s) = I m (A25) 


for s E <7, where N* and M* are as defined by equations (A16) and (A17), respectively. For a 
minimal (A, B, C, D) realization of G(s), a state- variable realization for an NRCF of G(s) is given 
by 


(A26) 


where S and R are given by equations (A 19) and (A20), respectively, 




A + BF 

BS” 1 /^ 

N 

M 

— 

C + DF 

DS -1 / 2 


F 

S " 1 / 2 


F = -S“ 1 (D t C + B t X) 


(A27) 


and X is the unique symmetric positive definite solution of the generalized control algebraic Riccati 
equation (GCARE) 

(A - BS _ 1 D T C) r X + X(A - BS _ 1 D r C) - XBS-^X + C T R _1 C = 0 (A28) 

Both the NLCF and NRCF realizations are minimal given a minimal (A, B, C, D) realization 
(Meyer 1988). 
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Appendix B 

Algorithms for the NLCF Robust Stabilization Problems 

In this section, computational procedures are presented for determining state- variable realizations 
of the compensators resulting from the optimal (7 = 7 min ) and suboptimal (7 > 7 min ) NLCF robust 
stabilization problems. 

Optimal Compensator 

From equation (23), construction of the optimal compensator involves the approximation of a 
function [— N, M] 6 RH 0 o by a completely unstable (all poles in the open complex right half-plane) 

ru v 

rational transfer function v .This type of problem is referred to as a Nehari extension problem. 

Its solution, as has been shown in section 6 of Glover (1984), can be constructed through the use of 
balanced realizations and Hankel norm approximations. In this paper, the version of Glover s theory 
found in section 6.6.2 of Maciejowski (1989) is employed. 

If G isapxm matrix, then [— N, Mj will bepx (p+m). Place an mx(p+m) null matrix beneath 
[— N, M] such that the resulting augmented square matrix is of size h = rn + p. A state-variable 
realization for the augmented matrix can be found through the use of equations (A 18) to (A22). 
The realization will have a system matrix with the same size as the system of a minimum realization 
of G, which is denoted by n. Let (A,B,C,D) be a balanced realization of the augmented matrix 
realization. The balanced realization will be such that 

AS + SA r + BB r = 0 (Bl) 


A ; S + EA + C 7 C = 0 


where 

£=[ CTlIr °1 (B3) 

.0 Si J 

Si - diag(rx r+ 1 ,<T r+ 2,---i CT Jj) ( B4 ) 

and r > 1 is the multiplicity of a\. In equations (B3) and (B4), the Hankel singular values of 
[-N, M] are ordered such that 


a i > a i+\ (i = 1 , 2 , . . . , n - 1 ) 


<71 =|| [-N,M] || H = || [N,M] || H 

Unique solutions to equations (Bl) and (B 2 ) exist because A is Hurwitz. 
Next, partition A,B, and C conformally with S to obtain 

An A l2 
A21 A22 



(B5) 


(B 6 ) 


and 

Define 



C = [GkC 2 ] 


r — S] — <7\lri—r 


(B7) 

(B 8 ) 

(B9) 
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and find U such that UU 7 = 1^ and 


(BIO) 


Bj = -CfU 


The matrix U may be found as follows. Let Cf have singular value decomposition UjEV^. 
From equations (Bl), (B2), and the construction of E, we have BiB, = Cj Ci, whereby UiEVj 
is a singular value decomposition of — Bj. Then, a solution of equation (BIO) is 


U = V1V2 (Bll) 

Thereafter, construct 

A = r _1 (<7iA 2 2 + SiA 22 Si -ctiC^UB^) (B12) 

B = r~ 1 (SiB 2 + aiC|’U) (B13) 

C = -C 2 E! - ctjUB^ (B14) 

D = — D + o\\J (B15) 


Glover (1984) shows that after deleting the lower m rows of C and D (corresponding to 

r u 1 * - - - - 

those initially added), a state-variable realization for L w \ in equation (23) is (A,B,C,D). 


Applying the definition in equation (A 16) indicates that the corresponding realization for 
is ( — A T , — C T , B^, D 7 ). In other words, 


U 

V 


= -B r (sl„_ r + A T )~ l C T + U T 


(B16) 


where 

A r € 

b t e /j(p +m ) x ( n ~ r ) 
C 7 e /j( n-r ) x P 

and 

D € 7?(P +m ) x P 


The matrices U and V are the upper rn x p and lower p x p portions of equation (B16), respectively. 
The optimal compensator K is given by equation (19) for which a state- variable realization can be 
constructed through application of equations (A4) and (A5). 

The system matrix corresponding to a state- variable realization of K will be on the order of 
( n — r), where r is the multiplicity of the largest singular value of [N,M]. Here, n is the order of a 
minimal realization of [N, M] and also the order of a minimal realization of G. Thus, the optimal 
compensator is of a lower order than G. However, even if G is strictly proper, the compensator will 
generally only be proper. (See eq. (B15).) If a strictly proper compensator is desired for a strictly 
proper G, the suboptimal compensator may be employed. 


Suboptimal Compensator 

Utilization of the compensator from the previous section always yields 


K 

I 


(I-GK)-'M - 1 Hoc = t ma X = T'min 


(B17) 
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If one wishes to specify that 7 > 7 min such that 


II 


K 

I 


(I - GK) - 1 M _1 


lloo < 7 


(BIS) 


a suboptimal version of the optimal NLCF robust stabilization problem can be posed as finding all 
compensators that simultaneously stabilize G and satisfy equation (B18). 

Let G have state-variable realization (A, B, C, D) and Z and X denote the solutions of 
equations (A 18) and (A28), respectively. Also let 


s = Im + D t D (B19) 

F=-S- 1 (B r X + D r C) (B20) 

and 

A r = A + BF (B21) 

Then, as shown in Glover and McFarlane (1989) or McFarlane and Glover (1990), all controllers 
satisfying the suboptimal robust stabilization problem for selected 

7 > 7min = (1- II [N, M] (B22) 

are given by 

K = (Ln$ + Li 2 )(L 2 i$ + L22) _1 (B23) 


where 


L,i 

! L12 


A c 

— 7 2 Wf 7 BS - 1 '' 2 

j 7 2 C _r ZC r R- 1 / 2 


I 

1 

= 

F 

s - 1 / 2 

j C _1 D r R“ 1/2 

L 21 

1 L 22_ 


C + DF 

DS -1 / 2 

| -C _1 R _1/2 


(B24) 


R = Ip + DD t 

C = (7 2 -l ) 1/2 
Wi = I n + (XZ - 7 2 I n ) 

and ^ is arbitrary in RH ^ so long as || # floo < 1 . 

The central (^ = 0), or lowest order, suboptimal compensator is 


(B25) 

(B26) 

(B27) 


K c = 

’ A c + 7 2 Wf T ZC T (C + DF) 

y 2 W~ T ZC T 


l 

to 

X 

-n T 


(B28) 


Comparing the orders of the optimal and central suboptimal controllers shows that the central 
controller always has an order equal to that of the original design model (G), whereas the order of 
an optimal compensator is n — r(r > 1 ). The central controller gives a strictly proper compensator 
for a strictly proper G, whereas the optimal compensator typically does not. 

McFarlane and Glover (1990) show that 


1 "hnin ~ ^max(ZX) 


(B29) 
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whereby Wi of equation (B 27 ) becomes singular as 7 -* 7min- As 7 -► 00, the central solution 
approaches 

[ A c - ZC r (C + DF) -ZC T 

B r X -D r 

which is easily recognized as a Linear-Quadratic-Gaussian (LQG) compensator construction. 
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Appendix C 
Numerical Data 

Matrices Defining Controller Design Model 

Numerical data for G = C(sli8 — A) _1 B + D for equation (49) are given in the following matrices 
for (A, B, C, D). 

A = 


Columns 1 through 6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 8.5437E-01 

0 

0 

0 

0 

0 

0 

— 8.7713E-01 

0 

0 

0 

0 

0 

0 

-9.5105E— 01 

0 

0 

0 

0 

0 

0 

— 2.1042E+01 

0 

0 

0 

0 

0 

0 

-2.2077E+01 

0 

0 

0 

0 

0 

0 

— 3.0155E+01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Columns 7 through 12 





0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 9.2432E— 03 

0 

0 

0 

0 

0 

0 

— 9.3655E— 03 

0 

0 

0 

0 

0 

0 

— 9.7522E-03 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-8.5712E+01 

0 

0 

0 

0 

0 

0 

— 1.1927E+02 

0 

0 

0 

0 

0 

0 

-1.3997E+02 

0 

0 

0 
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A = 


Columns 13 through 18 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 4.5872E-02 

0 

0 

0 

0 

0 

0 

-4.6986E-02 

0 

0 

0 

0 

0 

0 

-5.4914E-02 

0 

0 

0 

0 

0 

0 

— 9.2581E— 02 

0 

0 

0 

0 

0 

0 

— 1.0921E— 01 

0 

0 

0 

0 

0 

0 

— 1.1831E— 01 


26 


B = 


Columns 1 through 6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 9.0412E-02 

— 1.7503E-03 

-2.2032E-03 

— 7.2842E-01 

2.1216E— 04 

-7.4012E-01 

1.1157E+00 

— 2.6799E— 04 

6.2103E-02 

— 2.3142E— 04 

-1.1076E+00 

— 1.5399E-04 

5.9439E— 05 

— 7.5943E— 01 

1.9928E-04 

— 6.8777E— 02 

2.3805E-03 

— 2.5082E-01 

— 9.0841E— 01 

— 9.2293E— 04 

4.7562E-01 

2.9426E— 02 

1.2695E+00 

— 2.6966E— 02 

9.4147E— 01 

— 4.7622E— 02 

-9.0573E-01 

Columns 7 through 8 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


7.2043E— 01 

— 4.6856E-02 


2.8556E— 03 

— 6.2905E-01 


5.9586E-02 

5.9095E— 01 


— 2.7368E— 01 

1.0451E— 03 


3.9825E— 03 

— 1.1842E-03 


5.0834E— 04 

1.3185E+00 


— 7.0441E— 03 

1.1182E+00 


9.5012E— 01 

1.6953E— 02 


— 2.2581E— 02 

3.1281E— 01 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8.6983E— 04 
1.9792E— 04 
9.7330E— 05 
— 3.9490E— 02 
— 7.4550E— 01 
2.0583E-03 
4.0256E— 03 
— 1.0421E+00 
4.2744E— 02 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7.2204E-01 
2.8522E— 03 
5.9780E— 02 
-8.7477E-02 
— 3.0950E— 03 
8.2306E-04 
9.0757E-04 
— 3.8239E-01 
1.2593E-02 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7.6707E-02 
— 7.0804E— 01 
-8.8555E— 01 
2.3484E— 04 
-3.4298E— 04 
3.0302E-01 
— 4.3450E-01 
1.0193E— 02 
3.8921E— 01 
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c - 


Columns 1 through 6 


7.7245E— 02 

6.3892E-01 

-1.0611E+00 

1.4954E— 03 

— 1.8609E-04 

2.5487E— 04 

1.8823E— 03 

6.4918E— 01 

— 5.9063E-02 

— 7.4315E— 04 

— 1.7360E— 04 

— 9.2566E— 05 

-6.1689E-01 

-2.5017E— 03 

— 5.6854E— 02 

— 6.5536E— 02 

6.2104E— 01 

8.4221E— 01 

-6.1551E— 01 

-2.5047E— 03 

— 5.6670E-02 

4.0032E— 02 

5.5176E— 01 

— 5.6203E— 01 

Columns 7 through 12 


7.7862E+01 

-3.5096E+00 

— 1.3178E+02 

7.9107E— 02 

-1.5141E+02 

6.6658E+00 

— 4.0767E+01 

3.2162E+00 

1.2678E+02 

— 3.4504E-01 

1.2429E+02 

— 5.9830E+00 

— 7.7790E— 02 

4.5607E+01 

-1.7627E+00 

3.7242E+01 

-1.2157E+00 

-5.4479E+01 

6.0377E— 01 

— 1.1332E+02 

3.1607E+00 

— 9.5844E+01 

— 2.0220E+00 

— 4.3785E+01 

Columns 13 through 18 


1.0616E— 05 

— 2.7928E-06 

3.7768E— 03 

5.0808E— 02 

3.5683E— 02 

— 1.3072E— 04 

7.0638E-06 

— 9.3634E— 06 

1.3774E— 02 

1.8115E— 03 

3.5028E— 02 

— 1.1303E— 04 

4.0127E— 03 

1.4542E— 04 

-4.5198E— 05 

— 1.0773E-05 

1.6115E— 05 

-1.6640E— 02 

1.2554E— 02 

-1.8712E-04 

-2.7915E— 05 

— 4.7941E— 05 

5.5641E— 05 

— 7.2404E— 02 

D = 



Columns 1 through 6 


3.5007E+00 

-6.8878E-03 

— 6.5971E— 01 

— 6.8878E— 03 

3.4174E+00 

7.7127E— 03 

— 6.5971E— 01 

7.7127E-03 

1.6618E+00 

5.6295E— 03 

-7.1509E-01 

— 9.4997E— 03 

— 9.1940E— 04 

— 3.8808E— 01 

— 8.4551E— 04 

2.6140E— 01 

— 4.5194E— 03 

— 1.6658E— 01 

1.2391E— 02 

1.5061E+00 

— 8.6038E— 03 

3.1031E— 01 

8.2638E— 03 

4.1973E— 01 

Columns 7 through 8 


1.2391E— 02 

3.1031E-01 


1.5061E+00 

8.2638E— 03 


-8.6038E— 03 

4.1973E— 01 


— 9.8264E-01 

3.6531E— 03 


1.8407E— 01 

— 8.3024E— 04 


4.5186E— 03 

— 4.5921E— 02 


1.5008E+00 

1.2062E— 03 


1.2062E-03 

3.8341E+00 



4.8696E— 03 

-1.3122E— 03 

2.0740E+00 

2.3307E+01 

1.6766E+01 

-7.1785E— 02 

3.2403E— 03 

-4.3995E-03 

7.5636E+00 

8.3096E— 01 

1.6458E+01 

-6.2069E— 02 

1.8407E+00 

6.8328E-02 

-2.4820E-02 

— 4.9416E— 03 

7.5719E-03 

-9.1377E+00 

5.7589E+00 

— 8.7921E-02 

-1.5329E— 02 

— 2.1991E— 02 

2.6143E— 02 

-3.9760E+01 

8.3570E— 04 

6.8220E— 03 

-1.0881E— 02 

1.6178E— 05 

— 1.9870E— 06 

2.6135E— 06 

2.0365E— 05 

6.9316E— 03 

-6.0564E— 04 

— 8.0400E— 06 

-1.8536E— 06 

— 9.4918E— 07 

— 6.6740E— 03 

— 2.6712E-05 

— 5.8299E— 04 

— 7.0902E— 04 

6.6311E-03 

8.6361E-03 

— 6.6591E— 03 

— 2.6744E-05 

— 5.8109E— 04 

4.3310E— 04 

5.8914E— 03 

— 5.7631E-03 

8.4102E— 02 

— 3.2136E— 03 

— 1.1139E-01 

8.5446E— 05 

— 1.3864E— 01 

5.6342E-03 

— 4.4033E— 02 

2.9450E— 03 

1.0716E— 01 

— 3.7269E— 04 

1.1381E— 01 

— 5.0570E— 03 

— 8.4024E— 05 

4.1761E— 02 

— 1.4899E— 03 

4.0226E— 02 

-1.1132E— 03 

— 4.6047E-02 

6.5215E— 04 

-1.0376E— 01 

2.6716E-03 

— 1.0352E— 01 

— 1.8514E— 03 

— 3.7009E— 02 


5.6295E— 03 

— 9.1940E— 04 

2.6140E— 01 

— 7.1509E— 01 

— 3.8808E— 01 

— 4.5194E— 03 

— 9.4997E— 03 

— 8.4551E— 04 

— 1.6658E— 01 

1.6452E+00 

4.0543E— 01 

4.9757E— 03 

4.0543E— 01 

6.7897E— 01 

1.2671 E-03 

4.9757E— 03 

1.2671E— 03 

1.7236E+00 

— 9.8264E-01 

1.8407E— 01 

4.5186E— 03 

3.6531E— 03 

— 8.3024E— 04 

— 4.5921E— 02 
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Matrices Defining Truncated System 

Numerical data for AG = Cf(sl 32 - A*) _ 1 B* + D* for equation (50) are given in the following 
matrices for (A t , B t ,C t ,D t ). 

A t = 


Columns 1 through 6 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

— 2.0909E+02 0 

0 -3.1809E+02 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 0 
0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

— 6.3630E+02 0 

0 -6.4186E+02 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

o 0 

0 0 

0 0 

0 0 

o o 

0 0 

0 0 

o o 

-6.9828E+02 0 

0 — 7.6148E+02 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 
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Columns 7 through 12 

0 0 

0 0 

0 0 

0 0 

0 0 

o 0 

o 0 

0 o 

o 0 

0 o 

o 0 

0 0 

0 o 

o 0 

0 o 

o 0 

0 o 

o 0 

0 o 

0 o 

o 0 

0 o 

— 1.1948E+03 o 

0 -1.5075E+03 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 o 


0 0 

0 o 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

0 o 

0 o 

o 0 

0 o 

0 o 

0 o 

0 o 

0 0 

o 0 

0 0 

o 0 

0 o 

o 0 

o 0 

0 o 

— 1.5326E+03 o 

0 — 1.6531E+03 

0 0 

0 0 

0 0 

0 0 

0 0 

o 0 


o 0 

o 0 

0 o 

o 0 

o 0 

o 0 

o 0 

0 o 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

0 o 

o 0 

o 0 

o 0 

0 o 

0 o 

0 o 

— 1.7563E+03 o 

0 — 2.1457E+03 

0 0 

0 0 

0 0 

0 0 
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Columns 13 through 18 

0 0 

0 0 

0 0 

0 0 

0 o 

0 o 

0 o 

0 o 

o 0 

o 0 

o 0 

o 0 

o 0 

0 o 

0 o 

o 0 

o 0 

0 o 

o 0 

o 0 

0 o 

0 o 

o 0 

o 0 

o 0 

0 o 

0 0 

0 0 

— 2.7152E+03 0 

0 -3.1739E+03 

0 0 

0 0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 6.1546E+03 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 1.1213E+04 


1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.4460E-01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 1.7835E-01 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Columns 19 through 24 


0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 2.5225E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 2.5335E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 2.6425E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2.7595E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 3.4566E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 3.8827E— 01 
0 
0 
0 
0 
0 
0 
0 
0 
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Columns 25 through 30 

0 o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

l.OOOOE+OO 

0 

0 

0 

l.OOOOE+OO 

0 

0 

0 

l.OOOOE+OO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 3.9149E— 01 

0 

0 

0 

-4.0658E-01 

0 

0 

0 

-4.1908E-01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

l.OOOOE+OO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 4.6322E— 01 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000E+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 5.2108E— 01 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

l.OOOOE+OO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 5.6337E— 01 
0 
0 
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Columns 31 through 32 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

o 0 

o 0 

o 0 

o 0 

o 0 

o 0 

1.0000E+00 o 

o 1.0000E+00 

o 0 

0 0 

o 0 

o 0 

0 0 

o 0 

o 0 

o 0 

0 0 

o 0 

o 0 

o 0 

0 0 

0 0 

— 7.8451E — 01 o 

0 — 1.0589E+00 
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B t = 


Columns 1 through 6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 1.7116E— 01 

— 5.8049E— 03 

2.0869E-01 

— 2.9829E— 03 

— 3.5789E-01 

— 1.2875E— 03 

— 2.0542E— 01 

-7.2812E-02 

— 9.5226E-03 

— 8.3452E-01 

— 1.7241E— 02 

— 4.2512E— 02 

— 2.7311E— 01 

-1.0499E-01 

— 4.2532E— 02 

8.3909E-02 

-1.0430E+00 

-4.9171E-03 

3.0662E-01 

— 4.5815E— 02 

4.1271E— 01 

1.5476E— 02 

-5.8067E-01 

-3.6281E— 02 

1.1892E-01 

-3.9154E-02 

3.6029E-01 

— 2.0457E-02 

3.3652E-01 

1.9963E-01 

— 6.8681E-02 

1.0751E+00 

— 1.2244E— 01 

5.4524E-03 

5.4087E-01 

1.5533E— 01 

— 2.0571E-01 

2.7394E— 02 

— 9.8462E— 01 

-1.5396E+00 

-1.0585E-01 

— 1.7445E— 01 

— 1.0503E-02 

1.9177E-01 

3.6153E-02 

1.7533E— 01 

1.5125E— 01 

1.1805E— 01 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3.2349E— 03 

-1.0585E-03 

1.0457E— 01 

2.5253E— 01 

— 3.8650E-02 

— 3.3579E— 03 

7.2871E-03 

— 5.5209E— 02 

6.7861E— 02 

-2.1983E— 02 

-2.8225E-03 

2.8268E— 01 

— 5.3459E— 03 

— 5.4968E-02 

1.0531E— 01 

-6.5816E-02 

-4.8285E-01 

-3.5314E— 02 

— 1.2829E— 02 

6.0601E— 03 

-8.9994E— 02 

-3.9609E— 01 

— 1.8812E-01 

— 3.4468E— 02 

-9.2921E— 02 

2.3797E-02 

3.5179E-01 

2.5025E— 02 

-1.8725E-02 

-2.1458E-01 

— 1.4175E— 01 

6.5808E-02 

3.4202E-02 

1.0847E+00 

— 6.6816E-01 

2.2399E-02 

1.0100E-01 

— 4.6903E-02 

4.6728E-02 

-6.0191E-03 

— 1.5392E-03 

— 1.2358E-03 

6.2227E-01 

4.0364E-01 

5.9594E— 03 

7.5504E— 01 

— 2.2316E-02 

2.6324E-02 
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Columns 7 through 8 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

— 7.9081E— 03 
— 2.0767E— 01 
-4.1550E— 02 
— 5.3852E— 02 
— 6.0046E— 02 
— 4.0749E— 01 
1.0370E— 02 
1.6912E— 01 
— 1.0616E— 01 
-5.2492E— 01 
-2.2117E+00 
5.2024E— 01 
1.0856E— 01 
1.1115E— 01 
3.7390E— 01 
1.0209E— 01 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 2.1884E-01 
3.3597E-03 
1.1263E-01 
4.5806E-01 
1.5792E— 01 
— 4.2931E-02 
— 5.1029E-01 
8.8657E-02 
— 1.1353E+00 
-1.2112E+00 
3.0859E-01 
5.3702E-02 
— 3.1874E-01 
— 2.4079E-01 
1.3817E-02 
— 7.8110E-03 
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Ct = 


Columns 1 through 6 


3.5788E+01 

9.4882E-01 

1.3071E+02 

1.2138E+00 

1.1384E+02 

4.6330E+01 

-4.3635E+01 

4.0954E-01 

6.0592E+00 

— 6.7639E— 01 

-8.0327E+01 

— 4.6368E+00 

2.2132E— 01 

1.2294E+01 

3.5130E+ 01 

— 2.1865E+01 

1.0681E+00 

-4.3180E+01 

1.6535E+00 

6.6057E+01 

2.6438E+01 

4.5758E+01 

— 1.0687E+00 

— 7.1667E+01 

Columns 7 through 12 


— 3.6635E+02 

-2.3331E+01 

— 1.8226E+02 

5.4740E+01 

8.7538E+02 

6.0009E+01 

-4.9311E+02 

5.4695E+01 

— 5.5220E+02 

1.5328E+01 

5.9712E+02 

1.4241E+02 

— 7.2407E+00 

2.8360E+02 

— 3.6472E+01 

1.0753E+02 

5.1962E+01 

— 5.3917E+02 

-1.2390E+01 

— 2.5495E+02 

1.6271E+02 

6.0970E+02 

-1.3365E+02 

1.7400E+03 

Columns 13 through 18 


5.5855E+02 

4.8865E+03 

6.4641E+01 

— 7.4381E+01 

3.3595E+02 

— 1.1803E+03 

2.6735E+03 

5.5368E+02 

-2.2251E+02 

— 2.7424E+02 

1.9104E+01 

-3.8298E+03 

1.2735E+02 

4.8852E+00 

— 2.4842E+03 

— 1.2688E+02 

3.9223E+00 

-3.6677E+01 

— 2.9477E+02 

-3.5277E+02 

-2.3012E+03 

8.6546E+02 

7.6423E+02 

-8.5038E+01 

Columns 19 through 24 


5.1817E-02 

2.1143E— 01 

7.2169E-02 

1.8367E-02 

4.3680E-03 

2.7744E— 02 

2.4021E— 03 

1.0770E— 02 

1.1239E— 02 

-1.8382E— 03 

5.5694E-03 

1.4127E— 03 

1.3926E-02 

7.1508E-04 

1.4525E— 02 

— 1.7118E-02 

-7.1617E-02 

— 2.7828E-02 

1.0481E-02 

1.3643E-02 

1.5867E— 02 

— 2.8411E— 02 

— 1.1605E-01 

— 4.1730E— 02 

Columns 25 through 30 


— 4.6556E-02 

8.3174E-03 

2.8783E-02 

1.5328E-02 

— 1.3682E— 01 

— 4.5055E-01 

— 1.4105E— 01 

-8.1166E-02 

5.1312E-02 

3.6378E-02 

— 1.0175E— 02 

5.9405E— 02 

— 9.3163E-03 

7.6132E-03 

— 2.7579E— 02 

— 1.3772E-01 

8.7244E-02 

-1.4333E-02 

4.1561E-02 

2.1342E-01 

9.2688E— 01 

4.4446E— 01 

4.9245E— 01 

— 1.2932E-01 

Columns 31 through 32 


8.2397E-03 

-1.8566E-01 


-1.5045E-01 

-1.6016E-01 


— 2.8362E— 02 

-1.2500E-01 


— 4.8818E— 01 

— 7.9951E— 01 


-3.1666E-01 

2.3630E-02 


— 4.6752E-03 

-2.7874E-02 


— 2.9333E-01 

— 1.0810E-01 


— 1.0840E-02 

8.2711E-03 



5.3565E+02 
1.1066E+01 
2.7287E+01 
1.4110E+01 
1.8117E+00 
— 1.8144E+02 
3.4566E+01 
— 2.9401E+02 

3.3817E+01 

-5.5629E+02 

-3.3000E+02 

-4.1368E+01 

3.0954E+01 

3.5472E+02 

8.6773E+02 

2.0022E+03 

-1.9659E+03 
-1.6959E+03 
— 1.3237E+03 
-8.4660E+03 
2.5022E+02 
-2.9516E+02 
-1.1447E+03 
8.7582E+01 

-2.3155E— 02 
2.8782E— 01 
1.3569E— 03 
1.8162E— 02 
1.3324E— 01 
9.7449E— 03 
1.1245E— 01 
1.1847E-02 

— 2.5257E— 03 
— 2.5054E— 01 
— 7.1952E-02 
— 5.0245E— 01 
3.0951E— 01 
-1.0376E-02 
— 2.4099E-01 
— 2.4876E— 02 


1.9071E+02 
7.3312E+01 
2.9699E+01 
3.7329E+00 
3.8383E+01 
-7.3536E+01 
4.1929E+01 
— 1.1027E+02 

1.2062E+02 
— 1.8882E+03 
2.1504E+02 
2.4895E+02 
-1.1558E+02 
-6.0068E+01 
3.8844E+03 
— 5.4197E+02 

2.4750E-02 
8.3939E— 04 
-3.0177E-02 
— 4.6777E-04 
1.5306E— 04 

— 1.5121E-02 
1.1435E-03 
3.1644E-02 

-1.0599E-01 

1.5836E-02 

— 1.4266E-01 
4.4345E-03 

— 2.0947E-03 
3.1107E-02 
— 3.5845E— 03 
1.7639E-01 

1.0719E— 01 

— 1.4274E-02 
5.1307E— 01 

— 5.2629E— 02 
2.4440E— 02 
-2.4349E-02 
— 5.6568E-02 
1.6609E-01 


-6.3895E+01 
7.9423E+02 
3.7443E+00 
5.01 18E+01 
3.6768E+02 
2.6891E+01 
3.1030E+02 
3.2691E+01 

— 1.1699E+01 

— 1.1606E+03 
-3.3330E+02 
-2.3275E+03 

1.4337E+03 

-4.8062E+01 

— 1.1163E+03 

— 1.1523E+02 

5.3200E-04 
6.3830E-02 
2.2963E— 04 
— 4.5039E-02 
6.8932E-03 
5.9888E-04 
3.7038E-02 
— 5.9920E— 04 

— 6.0089E-03 
2.2546E-01 
1.4087E-02 
1.5379E— 01 
7.3041E— 02 
1.3383E-02 
— 6.5664E— 02 
— 3.4423E— 02 

8.6736E-01 
5.9633E-02 
9.8280E— 02 
3.3910E— 03 
8.6714E— 04 
6.9621E-04 
-6.2619E-02 
1.3565E-01 
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D t = 


Columns 1 through 6 
3.4067E+00 4.8934E-02 

4.8934E— 02 3.2066E+00 

6.7837E-01 3.4160E-02 

1.1984E-01 8.9253E— 01 

— 1.3647E— 02 4.1129E— 01 

— 2.8674E-01 1.0130E-02 

1.6819E-02 -1.7796E+00 

— 2.6614E— 01 -4.1046E— 05 

Columns 7 through 8 
1.6819E— 02 -2.6614E-01 

— 1.7796E+00 -4.1046E— 05 

1.1186E-01 — 6.1139E-01 

1.1030E+00 2.5951E— 02 

— 1.6368E— 01 — 1.3970E— 02 

4.2827E-03 3.2796E-02 

5.8696E+00 3.4138E-02 

3.4138E— 02 3.5793E+00 


6.7837E-01 1.1984E— 01 

3.4160E-02 8.9253E— 01 

1.4429E+00 1.8143E— 01 

1.8143E-01 2.3992E+00 

-3.4086E-02 -4.1076E-01 

9.6979E-03 2.0039E-02 

1.1186E-01 1.1030E+00 

— 6.1139E-01 2.5951E— 02 


— 1.3647E— 02 — 2.8674E-01 

4.1129E— 01 1.0130E— 02 

— 3.4086E— 02 9.6979E-03 

— 4.1076E— 01 2.0039E-02 

8.9345E-01 1.1979E— 02 

1.1979E-02 2.9147E-01 

— 1.6368E-01 4.2827E-03 

— 1.3970E-02 3.2796E-02 
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Table I. Open-Loop-System Characteristics of Phase 0 Model 



Eigenvalues 

Frequency, Hz 

n i ii7 

Mode 

Real 

ilmaginary 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

-4.622E-3 
-4.682E-3 
— 4.876E-3 
-2.294E-2 
— 2.349E-2 
-2.746E-2 
-4.629E-2 
— 5.460E— 2 
— 5.916E— 2 

— 7.230E— 2 
-8.918E-2 

— 1.261E-1 
— 1.267E— 1 
— 1.321E— 1 
-1.380E-1 
-1.728E-1 
-1.941E— 1 

— 1.958E— 1 

— 2.033E— 1 
— 2.095E— 1 
— 2.316E— 1 
— 2.605E— 1 
-2.817E-1 
-3.922E-1 
-5.294E— 1 

9.243E-1 

9.365E— 1 

9.752E— 1 

4.587E+0 

4.698E+0 

5.491E+0 

9.258E+0 

1.092E+1 

1.183E+1 

1.446E+1 

1.783E+1 

2.522E+1 

2.534E+1 

2.642E+1 

2.760E+1 

3.457E+1 

3.883E+1 

3.915E+1 

4.066E+1 

4.191E-i-l 

4.632E+1 

5.211E-1-1 

5.634E+1 

7.845E+1 

1.059E+2 

0.147 

.149 

.155 

.730 

.748 

.874 

1.474 

1.738 

1.883 

2.301 

2.838 

4.015 

4.032 

4.206 

4.392 

5.501 

6.180 

6.231 

6.471 

6.670 

7.372 

8.293 

8.966 

12.49 

16.85 
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Table II. Eigenvalues of K^(s) Optimal Compensator System Matrix With (A = f , 

■ * /n , max 

[(*•*- Q ) = (0-1,2, 0.5);e Amax = 0.6749] 


Real 

Complex j 

Real 

± Imaginary 

-1.2941 

-1.857E-1 

1.115E-1 

-.7197 

-2.983E-1 

7.433E-1 

-.7088 

— 1.342E — 1 

8.200E-1 

-.5482 

— 3.402E — 2 

4.579E+0 

-.5462 

— 2.898E-2 

4.685E+0 

—.5361 

-3.684E— 2 

5.468E-f 0 

— .5283 

-5.027E— 2 

9.245E+1 

— .5066 

— 6.042E— 2 

1.090E+1 

-.4936 

— 5.962E— 2 

1.183E+1 

— .4746 



-.4684 



-.4610 



-.4595 



-.3868 



-.3723 



1 1 


Table III. Closed-Loop Eigenvalues of Design Model Controlled by Optimal Compensator K(s) With M 

[(*, *, a) = (0. 1 , 2, 0.5); ex.max = 0.6749] 


Open-loop mode 

Complex ] 

Real 

± Imaginary 

1 

— 4.426E-2 

9.246E-1 

2 

-7.079E— 2 

9.485E-1 

3 

— 1.096E— 1 

9.645E— 1 

4 

— 2.716E— 2 

4.588E+0 

5 

— 2.658E— 2 

4.698E+0 

6 

— 3.237E-2 

5.491E+0 

7 

— 4.820E— 2 

9.258E+0 

8 

— 5.723E— 2 

1.092E+1 

9 

-5.968E-2 

1.183E+1 

10 

— 7.234E— 2 

1.446E+1 

11 

— 8.920E-2 

1.784E-f 1 

12 

— 1.26IE— 1 

2.523E-H1 

13 

— 1.268E— 1 

| 2.534E+1 

14 

— 1.321E— 1 

2.643E+1 

15 

— 1.380E-1 

2.760E+1 

16 

— 1.729E-1 

3.457E+1 

17 

— 1.941E— 1 

3.883E+1 

18 

— 1.958E— 1 

3.915E+1 

19 

— 2.033E— 1 

4.066E-H 

20 

— 2.097E— 1 

4.191E+1 

21 

— 2.316E— 1 

4.632E+1 

22 

— 2.606E— 1 

5.211E+1 

23 

— 2.818E— 1 

5.634E+1 

24 

— 3.923E-1 

7.845E+1 

25 

-5.295E— 1 

1.059E+2 


Damping ratio, < 


0.0478 

.0744 

.1291 

.0059 

.0057 

.0059 

.0052 

.0052 

.0054 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 

.0050 
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Table IV. Eigenvalues of K^(s) Suboptimal Compensator System Matrix With ex — 0.9f/t ,m&x 

[(fc,t,a) = (0.5, 1,0.1); ex, max = 0.6670] 


Real 

Complex 

Real 

^Imaginary 

a -0.1000 

-2.296E+0 

4.230E+0 

-.1612 

— 2.366E-1 

4.690E+0 

-.2945 

— 1.859E+0 

5.232E+0 

-.4769 

— 2.339E+0 

8.859E+0 

-1.8518 

-3.535E+0 

9.659E+0 

-3.1977 

-1.408E+0 

1.103E+1 

-6.8123 




“Repeated eigenvalues of multiplicity 8. 


Table V. Closed-Loop Eigenvalues of Design Model Controlled by Suboptimal Compensator K(s) With ex - 0.9fx.max 

[(At, i, a) = (0.5, 1.0.1); ex, max = 0.6670] 


Open-loop mode 


a l 

a 2 

3 

4 


Complex 


Real 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
If) 

17 

18 

19 

20 
21 
22 

23 

24 

25 


-2.610E-1 

-5.216E-1 

-1.767E-1 

or 

-9.510E-2 

— 4.760E- 1 
— 6.248E— 1 
— 9.416E-1 
— 4.921E-1 
— 8.062E-2 

— 1.056E— 1 
-1.266E-1 
-1.391E-1 
-1.336E-1 
-1.606E— 1 

— 1.747E-1 
-1.953E-1 
-2.038E- 1 
— 2.133E— 1 

— 2.122E- 1 
-2.334E— 1 
-2.623E-1 
-2.895E-1 
— 3.924E-1 
-5.295E-1 


^Imaginary 


Damping ratio, C 


8.872E — 1 
4.620E+0 

4.671E+0 

or 

4.653E+0 

5.480E+0 

9.237E+0 

1.086E+1 

1.178E+1 

1.447E+1 

1.785E+1 

2.522E+1 

2.538E+1 

2.643E+1 

2.766E+1 

3.458E+1 

3.884E+1 

3.919E+1 

4 .071 E-h 1 

4.194E+1 

4.634E+1 

5.212E+1 

5.639E+1 

7.845E+1 

1.059E+2 


0.2822 

.1122 

.0378 

or 

.0204 

.0865 

.0675 

.0838 

.0417 

.0056 

.0059 

.0050 

.0055 

.0051 

.0058 

.0050 

.0050 

.0052 

.0052 

.0051 

.0050 

.0050 

.0051 

.0050 

.0050 


"Not discernible from data. 
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Table VI. Eigenvalues of K, 4 (s) Suboptimal Compensator System Matrix With e A = 0 
[(k, i, a) = (0.5, l,0.1);f^ m(lx = 0.6670] 


Real 

Complex 1 

Real 

± Imaginary 

"-0.1000 

-7.687E-1 

5.570E-1 

-.3663 

-5.185E-1 

7.661E— 1 

-2.5134 

-1.039E+0 

4.555E+0 


- 1.568E-1 

4.661E+0 


— 9.240E— 1 

5.427E+0 


-1.188E+0 

9.153E+0 


-1.808E+0 

1.063E+1 

a i. . _ r 

-8.976E-1 

1.161E+1 


a Repeated eigenvalues of multiplicity 8. 


Table VII. Closed-Loop Eigenvalues of Design Model Controlled by Suboptimal Compensator 

K(s) With = 0 

[(A;, i, a) = (0.5, l,0.1);e Amax = 0.6670] 


Open-loop mode 

Complex 


Real 

± Imaginary 

Damping ratio, £ 

a l 

a 2 

3 

— 2.578E— 1 

8.945E— 1 

0.2769 

4 

— 5.448E— 1 

4.581E+0 

.1181 

5 

— 4.956E— 1 

4.647E+0 

.1060 


or 

or 

or 


-9.308E-2 

4.653E+0 

.0200 

6 

-4.497E— 1 

5.503E+0 

.0814 

7 

— 5.895E- 1 

9.275E+0 

.0634 

8 

— 8.884E- 1 

1.092E+1 

.0811 

9 

— 4.683E— 1 

l . 18 1 E+ i 

.0396 

10 

— 7.423E— 2 

1.446E+1 

.0051 

11 

-9.251E— 2 

1.784E+1 

.0052 

12 

— 1.262E— 1 

2.523E+1 

.0050 

13 

— 1.287E— 1 

2.535E+1 

.0051 

14 

- 1.324E-1 

2.643E+1 

.0050 

15 

— 1.417E-1 

2.762E+1 

.0051 

16 

— 1.731E— 1 

3.457E+1 

.0050 

17 

— 1.943E— 1 

3.883E+1 

.0050 

18 

-1.971E— 1 

3.916E+1 

.0050 

19 

— 2.048E— 1 

4.067E+1 

.0050 

20 

— 2.100E— 1 

4.192E+1 

.0050 

21 

-2.319E — 1 

4.633E+1 

.0050 

22 

— 2.608E— 1 

5.211E+1 

.0050 

23 

-2.828E- 1 

5.635E+1 

.0050 

24 

— 3.923E — 1 

7.845E+1 

.0050 

25 

— 5.295E— 1 

1.059E+2 

.0050 

a Not discernible from data. 
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(b) Equivalent compensator for unshaped plant (G). 
Figure 3. Loop-shaping procedure. 
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Singular value bounds 



Figure 4. Schematic of CSI Phase 0 Evolutionary Model. 
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Singular value bounds Singular value bounds 




(b) 16-mode truncated system (AG(s)). 

Figure 6. Unweighted open-loop nominal and truncated systems. 
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Frequency, rad/sec 

(a) Optimal compensator (K(s)) with e = e max . 



Frequency, rad/sec 

(b) Suboptimal compensator (K(s)) with e = 0.9e ma x- 
Figure 7. Unweighted optimal and suboptimal compensators (K(s)) with e max = 0.4417. 





Singular value bounds 3 Singular value bounds 











Figure 10. Weighted optimal compensator (K(s)) with a) 
K (*> = 


(0.1, 2, 0.5), e A = M, max = 0.6749, and 
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Frequency, rad/sec 

(b) Condition from inequality (40). 


Figure 11. Robustness conditions for inequalities (40) and (41) for weighted optimal compensator (K (s)) with 
* A ” f Amax ~ 0.6749 and K(s) = , ^ ( s ) • 
















Frequency, rad/sec 

(a) Condition from inequality (41). 



Frequency, rad/sec 

(b) Condition from inequality (40). 


Figure 15. Robustness conditions for inequalities (40) and (41) for weighted suboptimal compensator (K(s)) 
with £ A = 0-9O max, M.max = 0-667, and K(s) - 





Singular value bounds 


— i 

2J0-3 



1 0 ' 5 J-j 1 — i — 1 11111 1 > — * — » — * « i ■ i 1 , — i — I,,, 

10' 1 10 O i o 1 

Frequency, rad/sec 

Figure 16. Weighted suboptimal compensator (K A (s)) with (Ar,i,a) = (0.5, 1,0.1) and 



10-7 L 1 * — 1 — 1 — *- 1 1 * » 1 1 — i — i — i t > i ■ . ■ . . ■ 

10-1 100 10 1 , 

Frequency, rad/sec 

Figure 17. Weighted suboptimal compensator (K(s)) with = 0 and K (s) = 0 5 





Figure 18. Robustness conditions for inequalities (40) and (41) for weighted suboptimal compensator (K(s)) 
with e A = 0 and K(s) = ( s + o.i) K ^( s )- 
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2.5 



Figure 23. Input time history of actuator 7 for optimal compensator at ( k,i,a ) = (2, 2, 0.5). 



Figure 24. Output time history of accelerometer 7 for optimal compensator at (k,i,a) = (2,2, 0.5). 
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Figure 25. Input time history of actuator 8 for optimal compensator at (k,i,a) - (2, 2, 0.5). 



Figure 26. Output time history of accelerometer 8 for optimal compensator at (k,i,a) - (2, 2, 0.5). 
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Figure 29. Input time history of actuator 8 for suboptimal compensator at e A - 0.9c*,™* 
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Figure 30. Output time history of suboptimal accelerometer 8 at - 0.9r >4, max and 
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